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Abstract. We discuss metastable states in the mean- field version of the strong 
coupling BCS-model and study the evolution of a superconducting equilibrium state 
subjected to a dynamical semi-group with Lindblad generator in detailed balance w.r.t. 
another equilibrium state. The intermediate states are explicitly constructed and their 
stability properties are derived. The notion of metastability in this genuine quantum 
system, is expressed by means of energy-entropy balance inequalities and canonical 
coordinates of observables. 

PACS numbers: 05.30 ~ 74.20.F - 64.60.M - 42.50.L 
1. Introduction 

One of the most well-known phenomena in physics, is metastability. Phenomena like 
supercooled water or hysteresis in magnetic systems are easily observed and identified 
with metastability. But in spite of its clear appearance, the development of a complete 
theory of metastability is still a hard and unresolved problem, for an overview, see 
e.g. The basic concept of metastability could be formulated in the following 

way: a thermodynamic system is prepared in a special initial state different from the 
equilibrium state. If the initial conditions are suitably chosen, the system will not 
relax immediately to equilibrium, but it persists a longer period of time away from the 
equilibrium state, until some large fluctuation or an external disturbance occurs driving 
the system to equilibrium. The key problem is to formulate an expression for the size 
of the perturbation characterising the metastable regime 

Metastability has intensively been studied for classical models, and interesting 
results are obtained for the metastable relaxation in kinetic Ising models we will 

not try to give a full overview of this fleld but refer to one of these papers for more 
references. It was conjectured (see e.g. [Q) that the metastable evolution is governed by 
the growth of droplets of the stable phase in a background of the (initial) metastable 
phase. Small droplets are probable to disappear again, but bigger droplets tend to 
grow, ultimately driving the system to equilibrium. The lifetime of the metastable 
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phase should then be hnked to the probabihty of creating a droplet-excitation of critical 
volume 10. A rigorous result pointing in that direction was obtained in a paper by 
Schonmann and Shlosman 0] where they proved that the 'exit-time' could be expressed 
as a function of the equilibrium surface-tension of a Wulff-droplet of volume one. A 
Wulff-droplet is a typical droplet excitation, for which the shape is found minimizing 
the surface terms in the free-energy. This activity and success for classical models is in 
sharp contrast with the situation for quantum mechanical models where little is known. 
In a recent paper P , droplet states were constructed for the quantum mechanical XXZ- 
Heisenberg model. 

In this note, we want to develop ideas of @] in order to approach the phenomenon 
of metastability for quantum systems. We put forward that the underlying concept for 
exit-times and the metastable evolution is situated in the behaviour of non-extensive 
terms such as fluctuations of relevant observables, like the non-extensive (equilibrium) 
surface-tension of droplet excitations determines the metastable evolution in Ising 
models 0. For the BCS-model studied in this paper, we find a characterisation for 
metastability and define exit-times for different observables, expressed as a function 
of the equilibrium expectation values of the 'normal coordinates' of the observable 
under consideration. The metastable evolution at an arbitrary temperature between 
two extremal superconducting phases is studied. The evolution is driven by a semigroup 
with Lindblad generator P,^. In fact we consider a detailed balance dynamics between 
two different fixed phase states and realise the evolution from one equilibrium state to 
another equilibrium state with a different phase. The intermediate states are explicitly 
constructed and their thermodynamic properties are derived. These states are not 
invariant under the Hamiltonian evolution, but the correlation inequalities (section |^) 
are satisfied for the normal coordinates of observables. We present a simple criterion 
to distinguish observables exhibiting monotone relaxation or metastable relaxation, and 
in the latter case, an expression for the exit-time is given. The exit-time is the time 
necessary to leave the initial state on the basis of having reached the maximum value 
of the observable under consideration. In the BCS-model, the observables which are 



invariant under the gauge transformation group of the broken symmetry [jTO|-[T^, come 
over as relevant observables. They all exhibit metastable relaxation and all have the 
same exit-time. 

2. The BCS-model 

The strong-coupling BCS-model is described by the local Hamiltonians 

^ Af N 

where , and = a^a^ — a^a^ are the well-known Pauli matrices, if at acts on 
the Hilbert space (gj^^ Cf. 
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The equilibrium states studied here are elements of the set of (rf,/3)-KMS 
states in the thermodynamic limit (A^ ] oo) of this model, i.e. the states satisfying 
the equilibrium conditions of Kubo, Martin and Schwinger [13]: 

ujp{BA) = uop{AnpB) VA, BeB, (2) 

on the infinite tensor product algebra B = <S>iLi ^2 at an inverse temperature P and 
with the reversible Heisenberg dynamics: 

rJ.) = w- lim = e**^^. e"**-^'^. 
The extremal points of the set of equilibrium states are given by the symmetric product 



states 14 



uxi.) = lltTpx. (3) 

i=l 

on the infinite tensor product algebra B = <S)ili Pa is a 2 x 2 density matrix, given 
by the solutions of the gap-equation: 

^ exp[-(3hx] 

trexp[-phxy ^ ' 

where 

hx = ecr'^ - Act" - Act+, (5) 
and with order parameter A, given by, 

\ = tipx{cr+)=coxicJ+). (6) 
This equation can be transformed into the equation 

A(i_ltanh[/?A;]) = 0, (7) 



where k = ^/€^~+~\X\^■, {—k,k} is the spectrum of the effective Hamiltonians hx (H), 
which is independent of the phase of the order parameter A G C. It can readily be seen 
that this equation (|^ has always a solution A = 0. This corresponds to the normal 
phase state. Solutions with A 7^ exist if the following conditions are satisfied: 

/5^ = ^log(l±|). 

These solutions A 7^ 0, are understood to describe the superconducting phase states. 
The inverse temperature jS > (3c, fixes only |A|, the norm of the order parameter. The 
phase of the order parameter 0, defined by A = lAle*"^, (p G [0, 2tt], remains free to choose. 
This leads to an infinite degeneracy of the states in the superconductive regime and is 
due to spontaneous symmetry breaking |10|,|ll| . As this phase becomes important in the 




remainder of this article, it will from now on explicitly be labelled, i.e. we denote the 
order parameter (|^) as Ae*"^, with a phase (p ^ [0,27r], and norm A G M"*". Furthermore, 
we fix now a certain subcritical temperature (3 > (3^-, and hence the norm of the order 
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parameter, which will now be labelled by A G R"*". The superconducting pure phase 
states are therefore from now on denoted by uj^ instead of uj\x\ei't>- 

The mechanism of spontaneous symmetry breaking can be explicitly seen as follows: 
The Hamiltonian H^q is invariant under the norm-continuous gauge transformation 
automorphism group Q = {a^ \ip E [0, 2n]} on B, defined by the action: 

a^at) = e^V+ (9) 

On finite subsets A C N, these transformations are implemented by the unitaries 

lji> ^ g#QA/2 

where Qa is the (local) infinitesimal generator: 

Clearly a^^H^) = UfHjsiU"^ * = with {!,..., A^} C A, but this symmetry is broken 
in the equilibrium states of the superconducting phase, we have: 

^<t> [o^M)) = Ae^^*^'^^ ^ Ae^* = uo^). (10) 

In fact the gauge group Q establishes a relation between the superconducting pure phase 
states with different phase-factors: 

u^{a^{X)) = uj^+^{X), MX e B, and e [0, 27r]. 

In the following we learn that the observables X E B which are invariant under the gauge 
transformations, i.e. a^{X) = X, are relevant in order to determine the metastable 
evolution between two superconducting states. 

We conclude this paragraph by remarking that a common method to single out an 
equilibrium state with fixed phase 0, consists of adding a thermodynamic unimportant 
term to the local Hamiltonians (0), such as a vanishing external field: 

1 ^ 

-NY.^te-' + cr.e^'f, (11) 

i=l 

forcing the system in the limiting Gibbs state uj^ (D of the superconducting phase. This 
fixes the phase to 0. 



3. Thermodynamic stability 

An alternative way to characterise equilibrium states, is given by the correlation 
inequalities |TB[. These represent conditions, which have been proven to be equivalent 



to the P-KMS condition @, while their interpretation is related to the principle of 
minimum free energy. In other words, they are an expression for the thermodynamic 
stability of the KMS (equilibrium) states |jl3|: 
Energy-Entropy Balance Inequalities 

Let Tt{.) = lim^v e**''^^. e"**''^^ be the Heisenberg dynamics, (defined in a weak limit sense) 
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and (3 the inverse temperature. A state uo is a {rt, j3)-KMS state if and only if for all 
X e Dom{5) 

- tMX*6X) > u{X*X) log ^1^^, (12) 

with 6{.) the infinitesimal generator of the dynamics (rt). 

In this paper we concentrate on the BCS-model. Since this model is of mean- 
field type, we only need to consider product states (see section ||). For such states the 
correlation inequalities for operators X & B of the form X = Xi ^ X2 ® . . . ^ X^ 
with Xi,X2, . . . ,X„ G M2 follow from the correlation inequalities for Xi,X2, . . .X„. 
Therefore it is sufficient to consider only one-point operators, and (0) is reduced to: 
Energy-Entropy Balance Inequalities in the BCS-model 

A product state p on B = <S)^i ^2 is a {rt, [3)-KMS state for the BCS-model ^ if and 
only if for every one-point operator X G M2 the following inequality holds: 

(3p{X*[h,,X])>p{X*X)\ogP^^y (13) 

with hp the effective Hamiltonian in this state, i.e. hp = ecr^ — p{a~^)(J~ — p{a~)a^ . 



The interpretation of these inequalities is the following: the l.h.s. of (12) or 



(p!3|) refiects the change in energy if we alter the state with a 'perturbation X\ The 
r.h.s. of inequalities (|12|) or ( ]I3D originates from the corresponding change in entropy 
S{p) = — trplogp under the 'perturbation X\ For example, if the criterion ( [TB| ) is 
applied to unitaries U G M2, with U*U = 1 = UU*, and we substitute U for X in ([13|) , 
it is reduced to: 

piU*[hp,U])>0. 

This inequality expresses the fact that the the local state tr(p.) must have a lower 
internal energy than the perturbed states ti{UpU* .), while the entropy of these states 
remain unchanged: S{p) = S{U pU*). 

If a state a; satisfies condition ([l^) for an operator X, we say that u is stable under 
the 'perturbation X'. 

4. The perturbed states 



In classical kinetic models [0-^ the metastable evolution is introduced by a dynamical 
semigroup of dissipative maps (e.g. Glauber dynamics), satisfying the detailed balance 
conditions w.r.t. the asymptotic equilibrium measure. The notion of detailed balance 
and dissipative evolutions are generalised for quantum systems ||, Q . The quantum 
dynamical semigroups as well as the classical Glauber dynamics do share the same 
physical background, in the sense that they can be constructed as the result of a weak 
coupling of the system with a temperature reservoir system JTB|. 



A continuous 1-parametergroup {'^{t)\t > 0} of linear maps '-/{t) on the algebra 
B = <S)iLi ^2, is called a dynamical semigroup if for every t > 0, 7(t) is a completely 
positive, unity preserving map on B and 7(0) is the identity map. 
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Let L be the (densely defined) infinitesimal generator of such a dynamical 
semigroup, i.e. 

7(t) = e*^, Vt > 0. 

The dynamical semigroup is then called dissipative if the generator L is self-adjoint, i.e. 
L{A*) = L{A)* WA G B, and satisfies the following inequalities: 

L{A*A) > A*L{A) + L{A*)A, VA G B. (14) 

Let u he a. state on B, the dynamical semigroup (7(t))i>o is said to satisfy the detailed 
balance conditions w.r.t. the state lj if the following duality property holds: 

uj{XL{Y)) = uj{L{X)Y), VX, Y e B. 

Consider now one of the extremal superconducting phase states (^), we can construct 
a dynamical semigroup {7</,(t)|t > 0} with generator L^, satisfying the condition of 
detailed balance w.r.t. this state. Because we are dealing with product states, the 
generator of {'j(f){t)) is globally defined if it is defined on the local sites, i.e. 

L^:M2^ Ma. 

It is given by [|,§: 

L^{.) = ^exp[-/3(e. -6,)/2] ., i?.,] + [E*,, .]E,,) , (15) 

where Eij = stands for the matrix units in the base of eigenvectors 

G {—,+}} of h^, the effective Hamiltonian (^ of the (rf,/5)-KMS state uj^ @; 
iV'i) is the eigenvector corresponding to the eigenvalue ej of h^. It is immediately checked 
that is dissipative (0) and satisfies the detailed balance conditions in the state 
with fixed phase 0, i.e. 

u^iXL^Y)) = iu^{L4X)Y) VX, Y G M2. (16) 

The action of this generator (|15|) is naturally extended to operators X G i3 of the form 
X = Xi (g) X2 ® . . . ® X„ with Xi, X2, . . . , X„ G M2 by 

n 

L^{Xi O X2 ® . . . ® X„) = ^ Xi ® . . . ® L^{X,) O . . . ® X„. 

i=l 

Since uj(f, is a symmetric product state (0), the dissipativity and the detailed balance 
properties (0) are preserved for these more general operators and follow from the 
properties of their one-site factors. The dynamical semigroup (70(t))t>o is then defined 
by: 

70(t) = e*^^ Vt>0. (17) 

Remark that the detailed balance properties (|T6[) guarantee that u;^ is stationary under 
this dynamical semigroup, while any other (rt,/5)-KMS state with a different phase- 
factor is not invariant. 

Suppose now that our system is prepared in a /?-KMS state u^^j, at time t = 0. 
Based on ideas from classical models we apply at t > 0, an evolution (7</,(t)) (0) 
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with phase 0, and the system is forced to evolve accordingly. This implies that the 
system will ultimately relax to the equilibrium state cj^. The intermediate states ojt are 
introduced by: 

cut(.)=cu_<^(e*^^), (18) 

where the intermediate states ujt are constructed by applying the dynamical semigroup 
(7^(t))f>o (p^ with phase to the initial equilibrium state with phase —0. These 
states oJt are, by construction, again product states, and their density matrices can 
explicitly be calculated. They are understood to describe a metastable regime. For t 
small enough, the state Ut will still be close to the initial state oj^^, while if t becomes 
large enough, the system will be relaxing to the equilibrium state cu,^, where the words 
"close" and "almost" have to get a precise meaning. 

We proceed now with the explicit construction of the metastable states ut ([18|). 
The spectral decomposition of the generator (|TB|) and the dissipative evolution 70(t) 
is given in terms of the matrix units {Ef^) of the corresponding asymptotic effective 
Hamiltonian, which can be written as = kE'^^ — kEt^- We compute that: 

L<^(1)=0; L^{D) = -dD- L^E^^) = ~cE^^; L^{Et+) = -cEt+, (19) 

where d = 4 cosh(/3A;), c = 2 + 2cosh{f3k) and D = e'^^E'l_^ - e~^^'E^_. Hence, the 
expectation values for the operators ([19|) in the state Ut ([18|) are given by: 

^t(i) = 1; 

^,(Et)=e'*=^_^(4j; 
uj,{E^_^)=e-'^uj_^{Et^). 
From these equations (EUp and linearity, all expectation values in Ut, and hence the 



density matrix pt at time t, can be calculated. 
4.I. Exit-times and normal coordinates 

Remark that any observable X G M|" can be developed in its normal coordinates for 
the asymptotic equilibrium state (|T§), i-e. 

X = uj^X)! + a+(X) + a^iX) + aO(X). (21) 

The operators a'^{X) and ci^(X) are understood to be the creation, resp. annihilation 
operator of the normal modes determined by X, they are given in terms of X and 
the projection operators and Et^ on the eigenspaces of the asymptotic effective 
Hamiltonian h^, corresponding to positive, resp. negative energy: 

aiiX) = E^XEt.; (22) 

a-^{X) = Et^XEi^, (23) 

while a^{X) is the canonical constant of motion determined by X: 

aliX) = Et+XEt^ + Et_XEt, - u^{X)l. (24) 
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Indeed, aP^{X) is a constant of motion in the state tu^, since it commutes with the 
effective Hamiltonian and: 

^a;^(An(aJ(X))5) = iuj^{A[h^, ri(a° (X))]5) = 0, VA, 5, X G M^. 

Note that these operators |23|, all have expectation value zero in the asymptotic 
state ujff,. They are an expression for the fluctuations of X around its asymptotic 
equilibrium value uj^{X) (see also below). Using (|20|) and (^TJ) we get an expression 
for the expectation value of X in uot. 

UtiX) = ^^(X) +^_^(aJ(X) + a^(X))e-*^ + a;_^(a°(X))e-*'^. (25) 

The time limits t — > and t ^ oo of expression (|25| ) yield clearly the equilibrium 
expectation values uj_^{X) resp. uj^{X): 

\miu}t{X) =u;_<^(X); 

lim uJtiX) = uJipiX), 

t— >oo 

i.e. this time evolution describes the transition between the two states under 
consideration, i.e. the transition from cj.^ to u^. 

In general, the evolution of uJt{X) (|25| ) can express two types of behaviour as a 
function of the time t. This is easily derived from the analysis the functions 

)+ 



: t uJtiX), 



VX e Mf , 



fx it) : 

and their time derivatives. 

Let us illustrate this with two pictures (Fig. [4.1|) : 



(26) 





(a) fxAt)-t^MXi) 



(b) fxAt) ■■ t ^ CUt{X2) 



Figure 1. Typical pictures of monotone (a) and metastable (b) relaxation 



Firstly, uJt{X) can relax monotonically to its asymptotic value uj(f){X), and thus 



behaves qualitatively as shown in figure 4J.(a). This behaviour is met if the function 
fx{t) ( ^61) has only an extremum at t = 0. Calculating the time derivative of fxit) (|26|) , 
we see that this happens whenever 



c 

< -. 



^_^(a+(X))+^_^(a^(X)) - 



(27) 
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where 

c 1 + cosh(/3A;) 
1 ~ 2cosh{(3k) 

This condition is satisfied if e.g. the second and the third term on the r.h.s. of expression 
(p5|) have the same sign. The relaxation is exponentially fast and is determined by the 
constants c = 2 + 2 cosh(/5/c) and (1 = 4 cosh(/5A;) (|l^). 

The second possibility is that condition (^Tf) is violated, i.e. if 



^-<pWX)) ^ 1 + cosh(/3fc) 



cu_^(aJ(X)) + u^^{a-{X)) 2 cosh(/3fc) 

In this case the function fx{t) '■ t t-^ uJt{X) (H) has an extremum, reached at a time 
> 0, which may depend on the observable X. A typical picture of such a behaviour 
is shown if figure |4n](b). 

The interpretation of this extremum is the following. The time t^{X) is a time- 
scale which indicates where the expectation value uJt{X) leaves the metastable regime 
and the relaxation to the asymptotic equilibrium value uj^{X) starts. For t < t*(X), 
uJt{X) moves away from the original equilibrium value a;_(^(X), this behaviour with an 
increasing distance from equilibrium corresponds to the metastable regime. At time 
t = t*(X), this evolution comes to an end since uJt{X) is now at its maximum distance 
from the equilibrium value uj_^{X). For t > we see an other type of behaviour. 

The system is again approaching equilibrium, but uJt{X) is not returning to the initial 
equilibrium value lj-^{X), it relaxes now (monotonically) to its new equilibrium value 

Since the time t*(X) marks where uJt{X) leaves metastability and the relaxation to 
the asymptotic value is started, we call the exit-time for the observable X. The 

exit-time can be found calculating the extrema of function (|^), i.e. is the 

time t for which = 0, it is readily computed being: 




(29) 



In this expression, it is clear that the exit-time for the expectation value of an observable 
t^:{X) is determined by the constants c,d (^) and the ratio between the initial-state 
expectation values of the normal coordinates of X (^), (^) and (p^ . 

In general we can formulate the following statement: 
Metastable Relaxation for the BCS-model 

A system prepared in a initial product state uj-^ on B = ^'^i M2 will relax to a 
superconducting equilibrium state uj^jj (0j at subcritical temperature [3 (j^ according to 
the evolution induced by a dynamical semigroup ( ff^ with a dissipative generator 
(|7^, satisfying the quantum detailed balance conditions in state u^j,. The behaviour of 
this relaxation depends on the chosen observables X e M|".- 

• Monotone Relaxation 
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The relaxation of the expectation value of an observable X G Ml*^ is monotone if 
the following holds: 



< 



1 + cosh(/3fc) 



a7_<^(a+(X))+^_^(a-(X)) - 2cosh(/3A;) 

where aJ(X), a^(X) and a^^{X) are the normal coordinates of X in the equilibrium 
state uj^ as defined by equation (\2^, and k = Ve^ + A^, the positive eigenvalue of 
the effective Hamiltonians at inverse temperature (3 

Metastable Relaxation 

The expectation value of an observable X G M^" relaxes metastably if the following 
is true: 



> 



1 + cosh(/5fc) 



u_^{a+{X))+LU_^{a^{X)) 2cosh(/?A;) 

In this case we define the exit-time for the observable X , t*(X) as the time when the 
expectation value of X has past its extremal value and the relaxation to equilibrium 
starts. is then given by the expression: 



U{X) 



log 



2 cosh(/5fc) 



log 



2 cosh{/3k) - 2 \ 1 + cosh(/3A;) 

An interesting point of these considerations is that the decomposition (|2l|) and 
the construction of the exit-times is in terms of operators which are closely related to 
quantum fluctuations. Quantum fluctuations F^{X) are the limits n oo of operators 
Fn{X) defined by: 



1 " 

F4X) = —y^X,-uiX), 



one 



where X is a local observable and Xj is a copy of X, translated to site i. In |T7| , p^ 
has proved the existence of the limit: 

lim u {Fn{Xf) = Co {F^{Xf) , 

n^oo 

defining a dynamical system on the level of the algebra of fluctuations {F^|X}, and 
defining a state Cj on the algebra of fluctuations. Also a dynamics (ft) of fluctuations 
induced by the original one (r^), is defined by the formula 

nF^(x) =F^(nX). (30) 

In |T^,[T^ it is proved this dynamical system is always a bosonic system. This is worked 
out in great detail. We refer to |]T^|TB| for more details and more precise information 
on this subject. The point is that the exit-times (^) can formally be expressed in term 
of fluctuations, as follows: 



t*(X) 



d — c 



log - + log 

c 



cu_<^(F^(Xo)) 



Co_^{F^{Qx)) 



Metastability in the BCS-model 



11 



where 

1 " 

F^iQx) = lim 5^(aJ(X) + (X)), 



n . ^ 

n 



F^(Xo) = lim ^^(a5(X)) 



and is the hmiting state with phase —(f) on the fluctuation algebra. The 

thermodynamic limit in nominator and denominator is taken jointly. Although the 
mathematical formulation of this computation is far from completely coherent with the 



general theory of quantum fluctuations |T^,|T8[, we are tempted to conjecture that this 



construction might be the key to the understanding of metastability in a broader class 
of quantum systems. 

4-2. Invariant observables under the gauge transformation group 

Let us now continue the study of the metastable relaxation in the BCS-model between 
two equilibrium states, in particular the evolution from U-cf, to Ucf, for gauge invariant 
observables. As explained, these superconducting states are not invariant under the 
gauge-symmetry a^, ip G [0, 27r] (P) of the BCS-Hamiltonians (|I]), in particular: 

uj^ = o a2^. 

Observables X G M2 which are invariant under this gauge transformation group Q 
i.e. satisfying 

a^(X)=Xor [(T^X] = 0, (31) 
have the same expectation value in all superconducting phase states, hence 
uj^^{X) = uj^{X). 

Developing the observable X in normal coordinates w.r.t. cj,^ ( |2TD yields 

uj^^ (a+(X)) + uj_^ {a~iX)) + uj^^ (aJ(X)) = 0. 
Hence we can rewrite equation (|25|) yielding: 

uj,{X) = oo^iX) - uj^^ (aJ(X)) (e-*^ - e"*^) . 

This implies that the relaxation for observables invariant under the gauge transformation 
group Q (^ is of the metastable type. The transition between the metastable regime 
and the relaxation regime happens at the same moment i.e. the invariant observables 
a^{X) = X share the same exit-time (|29|) : 

U = i^^^ni^ = ^S^log ( i + . (32) 

d-c 2Vl-4A;2-2 ^ V2 / 

The time gives a time-scale for the transition between the metastable regime t < 
during which the system evolves away from equilibrium and the relaxation regime t > t^, 
where the system evolves towards the new equilibrium state. 
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4-3. Explicit computations 

Let us now illustrate the metastability in the BCS-model with a few special observables: 
using ( pO|) or (^) we calculate the expectation values in the intermediate states ([T8|) of 
some typical observables: 

ujt{a+e-''^ + a~e'^) = 2A (l - 2 sin2(0) (A^e"*'' + e'e"*'^) //t') ; (33) 
u}t{i(T^e~'^ - la-e'^) = 2A sin(20)e-*^ (34) 

In equations (^,0), the monotone exponential relaxation towards the state is clear. 
The limits t 0{oo) of the expectation values give again the equilibrium values in the 
states a;_<^(a;0), in particular one computes also: 

lim cut(a+) = Xe''^ = toJa^). 

Note that we have the following bound on the time evolution of the 'condensate': 
X(t) = \uJt{cr^)\ < X, indicating that the 'condensation' is suppressed in the intermediate 
states. 

The situation is different for the evolution of the expectation value of a^, the 
generator of the gauge group (P): 

uj,{an = -2e - 4sin2(0)^ (e"*^ - e''') ; (35) 

Clearly, this is an example of an observable invariant under the gauge transformation 
group Q (0), and its relaxation is metastable. In the time limits t — > 0(oo), uJt{<J^) (^) 
tends to the equilibrium value — 2e, but \u!t{<J^) — u±if,{a^)\ goes through a maximum, 
attained at t^, (^2]), the exit-time for invariant observables . 

4.4- Temperature Dependence 

To conclude this section about the relaxation behaviour, we give some remarks on the 
dependence on the temperature. At the critical point (H), i.e. if (T | Tc) or (/? | Pc), 
the order parameter A vanishes. This also implies the following equalities for T = T^. 

and the metastabihty effects disappears completely. The constants (c, d) , governing the 
speed of relaxation decrease to the values Cc = 2 + 2cosh(/?e) resp. dc = 4cosh(/5e). 
Since the expression for the exit-times t*(X) ([29| ) depends on the ratio of the expectation 
values uj_^{a^{X)) and co'_<^(a_j^(X)) + co'_0(a^(X)), the behaviour in limit (A 0) 
of t*(X) can be different for different observables X G M|". Both uj_^{a^{X)) and 
uj-(f,{a'^{X)) + uj-(f,{a^{X)) tend to zero at the critical point (A 0), but the speed 
at which they decay to zero can be different. This depend on the observable X under 
consideration. 
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If we consider the metastable evolution between two groundstates (T = 0,/3 = oo), 
we find that there are no intermediate states, since the constants c, d (p^ ) governing the 
relaxation speed blow up to infinity as j3 tends to infinity, hence 




uj^fj, for t = 0; 
uj^ for t > 0. 



This results from the fact that the dynamical semigroup which we consider (^7|) becomes 
trivial for groundstates. The detailed analysis of the dynamics when T — > indicates 
the existence of a non-trivial evolution only on a different time-scale. 

Hence, if one would like to compare the dynamics at different temperatures, one has 
to rescale the time with an appropriate, temperature dependent, scaling factor. This 
scaling and the resulting dynamics ask for an independent analysis and is out of the 
scope of our considerations here. 



5. Stability-instability properties of the intermediate states 

Here we consider the stability properties of the intermediate states ujt{-) ( |18|) on the 
basis of the correlation inequalities (|13D, i.e. we consider the equilibrium or stability 
conditions: 

Pu;t{X*[h,,X]) > u,{X*X) log i ^^l^xX*) ) ' ^^^^ 

where ht stands for the effective Hamiltonian in the state uJt- 

ht = ea' -uJt{(y+)<y- -uJt{<J-)<y+. (37) 

From section ^ we know that uot is a /5-KMS state if and only if equation (^) is 
satisfied for all operators X E B. In this section we analyse to what extend the 
intermediate states ([T8|) are still stable, i.e. we check for which operators X G M2 
these states uot satisfy the inequalities (^). The meaningful operators in this context 
are the matrix units in the base of the spectral decomposition of ht (|37[), denoted by 
{E(t)^^, E{t)^_, E(t)^+, E(t) }. Compute first the expectation values 

uJt{E{t)^^) = ^ + kt + eft/{2kt); ^^^^ 
cut(E(t)_+) = -Xtft/kt, 



with Xt = |a;t(cr"*")|, kt = a/c^ + A^, the positive eigenvalue of ht (0) and ft = 
4sin^(0)eA^ (e"*'^ — e"**^) /k"^. The expectation values of the other matrix units can 

be determined from the expressions (|3BD, since uJt{E(t)++) = 1 — uJt{E{t) ) and 

u;t{E{t)+.) = u;t{E{t)_+). 



5.1. uJt is not an invariant state 

A first observation is that the states uot, t > are not invariant under the Hamiltonian 
evolution: this is easily derived from the expressions (|38|), calculating 

ut{\ht, E{t)+.]) = 2ktUJt{E{t)^_) ^ 0. (39) 
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5.2. Stability 

The position Qx and momentum Px observables of a normal mode are constructed in 
the usual way, by means of the creation and annihilation operators: 

a+(X) = E(t)++XE(t)__; (40) 
aT{X) = E{t)__XE{t)++, (41) 

yielding the following expression for Qx and Px'- 

Qx = ^{at{X) + a~{X))- (42) 

Px =-^(a+(X)-a-(X)). (43) 

These observables {Qx,Px) satisfy the canonical dynamical equations: 

t[ht,Qx] =2ktPx; 
t[ht,Px] = -2hQx, 

justifying the name normal coordinates. For the sake of completeness we also give a°(X), 
the canonical constant of motion of X: 

a°(X) = E(t)++XE(t)++ + ^(t)__X^(t)__ - iUtiX)l (44) 

These definitions are analogous to the ones in equations (|2^), ( p3|) and, (|^) in section H 
with an important difference. In section ^, we need the normal modes for the asymptotic 
equilibrium state (|T6|), whereas the expressions (|40|), ( p]) and, (|4^ define the normal 
modes w.r.t. the intermediate state Ut- We can now formulate the following stability 
statements: 



The correlation inequalities ( pq ) are satisfied for all linear combinations aQx + bPx, 
with a,b & C and Qx and Px as in d^l j resp. ^4^). 
For such an operators the correlation inequalities yield: 

uJtiEit)^^)>uJt{E{t)++). 

Since kt + eft/ {2kt) > (|38D, this condition is satisfied for every intermediate state 

UJt. 

The constants of motion, e.g. a^{X) satisfy the correlation inequalities. 

An operator C G M2 is a constant of motion in uJt if it satisfies [ht, C] = 0, such 

operators can be written as C = aE{t)^^ + bE{t) with a, 6 G C. It is easy to 

check that the correlation inequalities (|36| ) are trivially satisfied. Moreover, this 
stability property holds for any operator C e M2 satisfying [hp, C] = in a general 
symmetric product state with density matrix p and effective Hamiltonian hp, see 
e.g. (H). 
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5.3. Instability 

The instability of the intermediate states is reflected in the following statements: The 
correlation inequalities are satisfied for the creation operators af{X) j^jdj), VX G M2, 
but not for the annihilation operators a^{X) 

Substituting the operators af{X) (|40|) resp. a^{X) (l4ll) for X in the energy-entropy 
balance inequahties (|36|) yields: 

if X = atiX), then ^^^^flH ^ (^^^ 

if X = ar(X), then > e^^^% (46) 

i.e. ci;t could only be stable under both ^^'"(X) and a^{X), VX G M2, if the equality holds. 
It follows from the gap-equation that these equalities holds only for equilibrium 
states. The intermediate states are stable under the creation operators (051), but not 



under the annihilation operators (|46|) . This can be seen as follows: using (|38D , it is 
readily computed that for all t > the strict inequality (|45| ) holds in the limit e — > 0. 
Suppose now that there exists a point in parameter-space (t', e'), with t' > 0, e' > 0, such 
that in that point the inequality (^) holds. By continuity of the state in parameter- 
space, we have continuity of the function / : (t, e) 1-^ ujt{E{t) ) — e^^'^*ti;t(-E(t)++). 

Since /(t', 0) < and /(t', e') > 0, there must exist a point where this function is zero, 
this would imply that the intermediate state at that point is an equilibrium state, which 
cannot be true, see e.g. (^). This proves that condition (^5]) is always satisfied in the 
intermediate states, while (|i6D yields that the stability is violated for the annihilation 
operators a~[{X). 

6. Outlook 

All results in this paper concern the BCS-model, and rely very much on the mean- 
field character 0] of this model: all the states under consideration are chosen within 



the set of symmetric product states, the extension of these results to other mean-field 



models ||12[ is straight-forward. However, this work should be considered as a prototype 
model of a scheme which can be generalised to bona fide interacting systems. The main 
argument for this is that metastability is formulated in terms of fluctuation operators 
and their dynamics. Relying on the general theory of quantum fluctuations [l^, |l^ 



these fluctuation systems are quasi-free or generalised free systems, and therefore we are 
confident that our results have a much wider validity far beyond the model considerations 
of above. We reserve this generalisation to interacting systems for a future contribution. 
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